Heralded state preparation in a superconducting qubit 
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We demonstrate high-fidelity, quantum non-demolition, single shot readout of a superconducting 
flux qubit in which the pointer state histograms can be resolved to below one part in 1000. In the 
weak excitation regime, continuous measurement permits the use of heralding to ensure initialization 
to a fiducial state, such as the ground state. This procedure boosts readout fidelity to 93.9% by 
suppressing errors due to spurious thermal population. Heralding also enables a simple, fast qubit 
reset protocol without changing the system parameters to induce Purcell relaxation. 



PACS numbers: 03.67.Lx, 42.50.Lc, 42.50.Pq, 85.25.-j 

Recent progress in superconducting qubits [T]-[4] has 
resulted in coherence times exceeding 10 /is [5[ |6] as well 
as the demonstration of multi-qubit algorithms [71 |8], 
further validating superconducting circuits as a viable 
platform for quantum information processing. The si- 
multaneous realization of a fast, high-fidelity, quantum 
non-demolition (QND) qubit readout, however, has thus 
far been difficult, with many schemes exhibiting either 
sub-unity visibility [9 , long measurement times [lOl [11] , 
or demolition of the quantum state [12\ Low power 
dispersive readouts based on circuit quantum electrody- 
namics (cQED) [13, 14 have demonstrated QND opera- 
tion [15] with near unity visibility, thus faithfully map- 
ping the qubit state to distinct frequency shifts of a mi- 
crowave resonator [16 . Traditionally, insufficient mea- 
surement sensitivity has limited the fidelity with which 
these pointer states can be resolved in this architec- 
ture. Recent advances in low- noise superconducting am- 
plifiers have provided a sufficiently high signal-to-noise 
ratio (SNR) [TTl [18] to resolve the state of the measure- 
ment cavity in a time much shorter than the qubit relax- 
ation time Ti, and thus achieve single-shot readout. 

In this Letter, we use a high speed readout based on 
a Josephson parametric amplifier (paramp) [19| to in- 
sert an additional measurement pulse before a qubit ma- 
nipulation and measurement sequence to verify that the 
quantum system is initialized in the ground state. Such 
heralding techniques are currently employed in other 
quantum information architectures, for example trapped 
ion [20, systems. With this technique, we effectively elim- 
inate state preparation errors due to the spurious excited 
state population commonly observed in superconducting 
qubits [21]. Furthermore, this method permits a rapid, 
deterministic reset of the qubit state — a particularly im- 
portant function in long-lived qubits where simply wait- 
ing for a time much longer than Ti is impractical. 

An overview of the qubit chip is shown in Fig. [l] (a). 
The flux qubit and quasi-lumped-element measurement 
resonator are fabricated from aluminum on a silicon sub- 
strate using double-angle shadow evaporation. The res- 



onator is formed by the parallel combination of an inter- 
digitated finger capacitor {Cr = 600 fF) and meander line 
inductor (L^ = 1.3 nH), providing a measured resonant 
frequency ujr/27r = 5.780 GHz. The resonator is cou- 
pled to the feed lines by planar coupling capacitors, 
which set the resonator linewidth f<i/27r = 9 MHz. A per- 
forated metal plane surrounds the resonator and forms 
the ground reference for an on-chip, co-planar waveguide 
(CPW) input for both the qubit manipulation and read- 
out pulses. 

The three-junction flux qubit (loop dimensions 3.8 x 
4.0 fim^) is inductively coupled to the resonator through 
a shared, 2.6-/im length of a 150-nm-wide constriction in 
the meander line inductor. The qubit is flux biased with 
a small superconducting coil mounted on the outside of 
the copper box containing the qubit chip. The cryopack- 
age is anchored to the 30 mK base temperature stage of 
a liquid-cryogen-free dilution refrigerator. A schematic 
of the cryogenic portion of the microwave measurement 
circuit is shown in Fig. [l](b). Qubit control and readout 
pulses are injected through the weakly coupled port of a 
directional coupler. The readout pulse is reflected from 
the resonator, acquiring a qubit-state-dependent phase 
shift. The reflected pulse passes through a series of cir- 
culators and a second directional coupler to the paramp, 
where it is amplified and reflected to the output port. 
The readout signal is further amplified by cryogenic and 
room-temperature amplifiers before it is down-converted 
to zero frequency and digitized in 10-ns increments. 

The paramp is flux biased to match the qubit read- 
out frequency and operated in phase-sensitive mode, 
wherein a strong pump tone is applied at the same fre- 
quency as the readout pulse. In this operating mode 
the phase of the reflected pump depends very sensitively 
on its power. If a small input signal is added in phase 
with the pump signal, the resulting small change in net 
power leads to a large phase shift. This form of am- 
plification theoretically adds no additional noise to the 
measurement, allowing the overall noise temperature of 
the system to remain close to the standard quantum 
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FIG. 1. Experimental setup. (Color online) (a) False color 
SEM image of the qubit and readout resonator with a mag- 
nified inset of the qubit. (b) Path of the readout signal at 30 
mK. 



limit Tq = huj/2kB = 139 mK. Previous research has 
demonstrated noise temperatures close to this limit [19]. 
All measurements in the current work were made with 
the paramp providing 20-dB gain with an instantaneous 
bandwidth of 7 MHz. 

We use spectroscopy to determine the flux qubit pa- 
rameters A/h = 6.15 GHz (energy splitting at the de- 
generacy point) and Iq = 204 nA (circulating super- 
current in the qubit loop) [^. The coupling strength 
between the qubit and resonator is g/27r = 105 MHz. 
In this geometry, the effective coupling strength is given 
by g' = ^sin(^), where sin(6>) = A/hujQi] co'oi/27r is the 
qubit transition frequency. We extract the total qubit- 
resonator mutual inductance Mt = g/Iqlr = 6.3 pH, 
where = ^/huJr/Lr = 54 nA is the rms zero-point 
current in the resonator [23 . The qubit is flux biased 
away from the degeneracy point to cJoi/27r = 7.80 GHz, 
corresponding to a detuning S/27r = 2.02 GHz. At this 
bias the dispersive shift is 2x/27r = 7.4 MHz, producing 
a 150° phase shift between the two qubit pointer states. 
We measure a relaxation time Ti ^ 1.8 /is, achieved at 
the expense of a reduced T| = 55 ns away from the flux- 
noise-protected degeneracy point. 

We first determine the measurement SNR, in particu- 



lar the degree to which we can separate the pointer state 
histograms. Our measurement protocol involves energiz- 
ing the readout cavity with an average photon population 
n. After 90 ns, when the cavity has essentially reached 
its steady state, we digitize the homodyne signal at 10^ 
samples/s. Figure |2] (a) shows a portion of the readout 
signal as a function of time for three separate measure- 
ment traces at n = 14.6 photons. The dashed line indi- 
cates the discrimination threshold voltage for associating 
the cavity response with the ground or excited state of 
the qubit. The system SNR is sufficiently high to enable 
single-shot discrimination of the readout pointer states in 
a single 10-ns integration bin. Quantum jumps from the 
excited to ground state of the qubit are observed as an 
abrupt downward jump in the signal. The characteristic 
decay time extracted from an exponential fit to the distri- 
bution of downward jump times agrees with the value of 
Ti obtained from the ensemble averaged measurements. 

In Fig. [2] (b), we present histograms as a log-linear 
plot of many ground and excited state preparations cal- 
culated immediately after the cavity has come to equi- 
librium. The presence of the smaller, secondary peak in 
each distribution indicates partial contamination of the 
pure states and leads to fidelity loss. The histograms 
are well separated but there is a small number of er- 
rant counts, that is, counts at values of the homodyne 
voltage not centered about the bimodal peaks. We at- 
tribute these counts to instances where a transition from 
one qubit state to the other occurs during the measure- 
ment. These transitions can be induced by noise or by 
Ti relaxation. 

We use these state histograms to calculate the read- 
out fidelity F = 1 — Pq — Pi where Pq and Pi are the 
fraction of error counts observed for a fixed discrimi- 
nation threshold voltage. The uncorrected readout fi- 
delity as a function of readout power is presented in Fig. 
|3] The '10-ns' fidelity series is calculated from the his- 
tograms of the homodyne voltage values from a single 
10-ns digitization bin, as described above. The 'inte- 
grated' data series is calculated by integrating the entire 
readout signal weighted with an exponentially decaying 
filter [24]. The optimal filter time constant is empirically 
determined for each readout power to maximize the fi- 
delity. At very low power, the '10-ns' fidelity is severely 
suppressed because the system SNR is low. A substantial 
fraction of the lost fidelity can be recovered through in- 
tegration. As expected, fidelity increases with increasing 
readout power and the two methods produce identical re- 
sults above n ~ 10, where the '10-ns' readout histograms 
become well-separated. The highest measured fidelity is 
91.1 ± 0.4% at n = 37.8. In Fig. |3]we also plot the qubit 
relaxation time Ti of the ensemble averaged signal dur- 
ing readout after an excited state preparation. As the 
power is further increased above n = 100, the fidelity 
and Ti both decrease as a result of increasing readout 
backaction. Below n = 14.6, Ti during readout always 
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FIG. 2. Qubit readout at n = 14.6 photons. (Color online) 
(a) Three individual qubit readouts starting at t = 0: excited 
state (solid red line), ground state (long blue dashes), abrupt 
quantum jump from the excited to ground state (short green 
dashes), (b) Log- linear histograms of many ground and ex- 
cited state preparations, with each histogram bin count nor- 
malized to the total number of counts. 
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FIG. 3. Raw fidelity and Ti during readout vs. readout power. 
(Color online) Fidelity curves measured using a single 10-ns 
acquisition (blue x) and using integration (red +) merge to- 
gether when the readout pointer state histograms completely 
separate. The relaxation time Ti measured during continuous 
readout monitoring (green o) is referred to the right axis. 



exceeds 1.5 /is; the structure in the data in this regime is 
likely due to frequency-dependent variation in the elec- 
tromagnetic environment seen by the qubit as the Larmor 
frequency is Stark shifted with increasing readout power. 

We now examine the measurement histograms in detail 
to investigate sources of fidelity loss. In particular, we as- 
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FIG. 4. Pure readout pointer state distributions. (Color on- 
line) A two-point correlation procedure is used to exclude 
qubit transitions and create pure ground and excited state 
distributions, represented by log- linear histograms, with each 
histogram bin count normalized to the total number of counts. 
By totaling the counts on the wrong side of the voltage dis- 
crimination threshold, the distributions are determined to be 
separated to better than one part in 1000. 



sess the degree to which the pointer state histograms can 
be separated. In Fig. [2](b), we postulated that counts 
found between the two readout pointer state peaks were 
due to quantum jumps during measurement. To discard 
these events and create pure- state distributions, we im- 
plement a two-point correlation procedure. The system 
is probed at times tA and t^, as shown in Fig. [2] (a). 
We retain measurements only when both of these read- 
ings return the same value for the qubit state. Thus, 
the readout exhibiting the abrupt quantum jump is ex- 
cluded. The 160-ns time difference between these points 
is several times longer than the response time set by 
the ~ 10 MHz system bandwidth, thus ensuring mini- 
mal autocorrelation between the signals at Ia and t^, 
and nearly random distributions dit t = tn (subscript H 
denotes histogram). The resulting histograms dX t = tn 
with these outliers removed are shown in Fie;. Hfor ^ 10 
ground and excited state readouts. Using the discrimina- 
tion threshold shown in the figure, we observe only 108 
false counts, thus allowing the pointer state histograms 
to be separated to one part in 1000. Consequently, the 
finite measurement SNR does not contribute to the ob- 
served fidelity loss. This procedure, however, does not 
correct for those rare events in which the system jumps 
twice. Unambiguously distinguishing these events is dif- 
ficult due to the limited system bandwidth, which could 
be increased in future experiments. 

To investigate other sources of fidelity loss, we herald 
the ground state by inserting a fast measurement pulse 
prior to qubit manipulation, as shown in the pulse se- 
quence of Fig. [5] (a). With the qubit in thermal equilib- 
rium with the environment, we energize the readout and 
extract the qubit state dX t = ts (subscript S denotes 
selection). If the qubit is determined to be in the ground 
state, no correction is applied; otherwise the subsequent 
readout is discarded from the total record. Events where 
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FIG. 5. Heralded state preparation at n = 14.6. (Color on- 
line) (a) Pulse sequence used for heralded state preparation. 
A pure ground state is heralded at t = ts. (b) Log- linear raw 
readout histograms for ^10^ excited and ground state read- 
outs (solid lines) compared with the histograms generated af- 
ter heralded state preparation (dashed lines). The histogram 
bin counts are normalized to the total number of counts. An 
increase in fidelity from 91.0% to 93.9% is reflected in the de- 
creasing size of the smaller peak of the bimodal distributions 
which are responsible for fidelity loss. 



the qubit is spontaneously found in the excited state can 
be attributed to either remnant thermal population or 
readout-induced excitation of the ground state [25] dur- 
ing the heralding pulse. With this procedure in place, 
we again prepare ground and excited state histograms, 
and observe a fidelity improvement of 2.9 ± 0.2%. We 
compare the log-linear readout histograms both with (se- 
lected) and without (raw) heralded state preparation in 
Fig. |5](b). The reduced overlap of the histograms with 
heralded state preparation is demonstrated by the de- 
creased size of the smaller of the bimodal peaks in both 
the excited and ground state distributions. This proce- 
dure corrects for the remnant thermal population of the 
excited state. 

To separate out the contributions of the remaining 
sources of fidelity loss, we examine long time traces 
(^ Ti) with the system prepared in the ground state. 
We record and analyze the statistics of individual tran- 
sitions between the qubit states, extracting the average 
readout induced excitation (F^) and decay (F^) rates. 
The calculated fidelity loss during measurement due to 
is 0.2%, while the contribution of F^ is negligible. 
The remaining, unaccounted-for sources of fidelity loss 
are estimated to contribute 1.5%. We speculate that this 
loss is explained by errors in excited state preparation, 
namely 7r-pulse imperfections and drifts of the Larmor 



frequency due to local magnetic fiux variations. In Table 
|lj all of the sources of loss and the method by which they 
are calculated are listed. We note that the dominant loss 
mechanism is energy relaxation, and calculate its contri- 
bution by comparing the readout time and the measured 
value of Ti . Currently, there are flux qubits with demon- 
strated relaxation times on the order of 10 /is [5l [26] : 
in our system, a similar Ti would reduce the associated 
fidelity loss to 0.7%. Moreover, one could potentially 
employ cavities engineered to suppress the spontaneous 
decay of the qubit [27 . This level of coherence, coupled 
with heralded ground state preparation, high-precision tt 
pulses [28 and stable magnetic flux, should readily en- 
able readout fidelities in excess of 98% within this low 
power dispersive measurement architecture. 



TABLE I. A budget of fidelity loss at n = 14.6. The measured 
fidelity is 91.0 ±0.4%. 

Source of loss Fidelity loss (%) Calculation method 



Ti decay 4.4 ± 0.3 

thermal population 2.9 ± 0.2 

Ft 0.2 ±0.1 

F^ 0.0 ±0.2 
SNR < 0.1 

remaining 1.5 



measured Ti 
heralding 
individual jumps 
individual jumps 
pointer state overlap 



Finally, we discuss a fast qubit initialization procedure 
based on heralding of the ground state. If the excited 
state is measured during the probe pulse, the state can 
be immediately reverted with a control pulse and rapidly 
checked again. This procedure eliminates the need to 
change the detuning 5 ^7] to induce Purcell relaxation. 
The ideal fidelity of this active reset to the ground state 
is ultimately limited by how the relaxation time of the 
qubit compares with and is given approximately by the 
simple formula Freset = e~^^l^^^ = 93% for our experi- 
mental parameters. We measure a reset fidelity of 82%, 
which is close to the reset fidelity we expect to achieve 
(84%) after correcting for the fidelity loss ascribed to ex- 
cited state preparation errors (1.5%) and an extra 140 ns 
of delay present in the pulse sequence used in this par- 
ticular experiment (8.1%). The reset procedure could be 
iterated to achieve enhanced reset fidelity, but for qubits 
with relaxation times greater than 10 /is, a single itera- 
tion of this procedure would achieve reset with close to 
99% fidelity, assuming that the leftover fidelity loss we 
attribute to imperfect state preparation could be elimi- 
nated. 

In conclusion, we have demonstrated a fast, analog 
cQED readout to measure the state of a flux qubit with 
high single shot fidelity. We use the QND nature of 
this readout at low powers to demonstrate heralded state 
preparation. This procedure allows us to obtain a de- 
tailed account of the sources of fidelity loss in the read- 
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out protocol. Moreover, errors due to imperfect ground 
state initialization can be readily identified and effec- 
tively eliminated through post-selection. Extending this 
further, these errors could be actively corrected in real 
time if a TT pulse were to be triggered based on the out- 
come of the heralding readout. The dominant remaining 
source of fidelity loss in our experiment is due to the short 
energy relaxation time of the qubit. For Ti > 10 /is with 
perfect state preparation, we predict a fidelity of greater 
than 98% in the current readout architecture. We also 
note that this simple heralding procedure immediately 
allows for rapid reset of the qubit without the need to 
change bias parameters. Our readout architecture pro- 
vides a new tool to study readout-induced backaction in a 
qubit with a high degree of anharmonicity, complement- 
ing previous work on transmon qubits [25l |29] . 
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